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Abstract. Recently, the modified BFW model on random graph [Phys. Rev. Lett., 106, 115701 (2011)], 
which shows a strongly discontinuous percolation transition with multiple giant components, has attracted 
much attention from physicists, statisticians and materials scientists. In this paper, by establishing the 

J: 

theoretical expression of evolution equations on the modified BFW model, the steady-state and evolution 
process are analyzed and a close correspondence is built between the values of parameter a and the number 
of giant components in steady-states, which fits very well with the numerical simulations. In fact, with 
the value of a decreasing to 0.25, the error between theoretical and numerical results is smaller than 4% 
and trends to rapidly. Furthermore, the sizes of giant components for different evolution strategies can 
also be obtained by solving some constraints derived from the evolution equations. The analysis of the 
steady-state and evolution process is of great help to explain why the percolation of modified BFW model 
is explosive and how explosive it is. 
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1 Introduction flow through a disordered porous medium. The research 

in percolation is not only of academic interest but also of 

Percolation is a classical model in statistical physics, prob- considerable practical value. During the last five decades, 

ability theory, materials science, complex networks and percolation theory has found a broad range of applica- 

epidemiology, which is initiated as a mathematical frame- tion in epidemic spreading, porous media, robustness of 

work for the study of random physical processes such as networks to attacks, etc [IJ[2j|3]. 
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Percolation has been studied on various topological Furthermore, K. J. Schrenk et al. |29j generalize the re- 
structures such as scale-free network, lattices with differ- suits to the lattice with different dimensions, 
ent dimensional, random graph, etc. Taking percolation So far, although one-dimensional and mean-field per- 
model on the Erdos-Renyi random graph (ER model) as eolations have been solved theoretically, the others still 
an example, this model is one of the most simple and clas- remain on researching, especially the explosive percola- 
sical models that undergo a phase transition of a emerging tion, which attracts much attention from physicists, statis- 
giant component. Typically, percolation phase transition ticians and materials scientists. As a typical member of 
is considered as a robust second-order transition until a explosive percolation, the percolation threshold of BFW 
recent work by D. Achlioptas, R. M. D'Souza, and J. model has been analyzed, but many other properties are 
Spencer [3], in which they propose that the phase tran- still not clear, which drives us to investigate the evolution 
sition of some certain Achlioptas process is discontinuous process of BFW model with both simulation and theoret- 
and call it explosive percolation. This interest phenomenon ical method. 

leads to intensive studies on the other models like scale- This paper is organized as follows. In section 2, we in- 
free network [5]|6|, local cluster aggregation model [Jj and troduce the BFW algorithm with parameter a in details; 
lattices [SUH] • More recently, it has been demonstrated that the mathematical expression of BFW model is established 
all Achlioptas processes have continuous phase transitions and analyzed in theory, by which, we obtain the steady- 
in the mean- field limit [TUl[mn2lll3U14[[T3] . But some state condition and evolution regulation of BFW model for 
other kinds of models, which have different and special any a. In section 3, by analyzing the evolution procedure 
rule of evolution, such as triangle rule [TB] , largest cluster of BFW model, we find the relationship between param- 
rule [17], etc [T51ll9U20U21l[2^1l23U24l[23] . have been ana- eter a and the steady-state, that is when a € (^j, 
lyzed in details and indeed exhibit explosive percolation. BFW algorithm must stabilize with m giant components, 

for any m 6 N + . Furthermore, size of these components 



must satisfy some constraint equations which are given in 
our paper. 



In particular, the BFW model on random graphs, orig- 
inally introduced by Bohman, Frieze, and Wormald [26 , 
is similar to Achlioptas processes but more restricted. The 
recent work of W. Chen and R. M. D'Souza [2"71l2"g] . shows 
a strongly explosive percolation with multiple giant com- The BFW model on random graph is firstly introduced 



2 Dynamical behaviors of BFW(a) model 



ponents in BFW model. It is also shown that with smaller by T. Bohman, A. Frieze, and N. C. Wormald [2B], aim- 
parameter values, the transition will become more explo- ing to choose a subset A € {ei,e2, ■■■,e,2t\ with \A\ = t 
sive and the number of giant components will increase, such that for t as large as possible the size of the largest 
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component in G — (n, A) is o(n) (i.e. G does not contain 
a giant component); here n denotes the number of nodes 
and {ei , e%, ...} are the sequence of edges chosen uniformly 
at random from the edge set of complete graph; A repre- 
sents the set of accepted edges (initialized to A = 0) and 
t = \A\ represents the number of accepted edges. 

According to the BFW model, one of the sampled 
edges is considered at each step, and either accepted to the 
graph or rejected provided that the fraction of accepted 
edges is never smaller than the decreasing function g(k), 
which is asymptotically approaching the value 1/2. If tak- 
ing u as the total number of sampled edges, the fraction of 
accepted edges is represented by t/u; k denotes the stage 
and the function g(k) = l/2 + l/v / 2fc- This model is much 
similar to Achlioptas process and shows that the thresh- 
old of a giant component is t — c*n where c* satisfies a 
certain transcendental equation and c* 6 [0.9792,0.9793]. 
This result has been verified by theoretical methods [35] 
and simulations [27] . 

Recently, the BFW model is extended to BFW(a) one 
and analyzed by W. Chen and R. M. D'Souza [27J with 
modifying the function g(k) to a + 1/ \/2k. It is shown that 
multiple giant components appear in a strongly explosive 
percolation transition. Furthermore, with the value of a 
decreasing, the threshold will delay and the phase tran- 
sition will be more explosive. In the following sections, 
we will discuss the reason and properties of these phe- 
nomenons and provide theoretical analysis. 

For the theoretical analysis on the BFW(a) model, the 
BFW (a) algorithm is shown as follows: 



algorithm BFW(a) 

1 begin 

2 A = 0; 

3 k = 2; 

4 t = u = l; 

5 e is a randomly sampled edge; 

6 while(£ < 2n) 

7 begin 

8 1= Maximum size of component in A[J{e}; 

9 if (/ < k) 

10 A = A(J{e}; 

11 t = t + l; 

12 u = u + l; 

13 sample an edge e randomly; 

14 else ii(t/u < a + l/y/2k) 

15 k = k + l; 

16 else 

17 u = u + 1 ; 

18 sample an edge e randomly; 

19 end 

20 end 

2.1 Evolution analysis on the BFW(o) model 

To analyze the evolution of BFW(q) model and its steady- 
states, we consider the following variables: k, t, u and n 
possess the same meaning as they are in the BFW(a) al- 
gorithm; m represents the number of components; Ci de- 
notes the fraction of ith largest component. 
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In BFW(a) algorithm, there are three cases when an (i t + 1) and the components Cj and Cj merge together 
edge is sampled: (m <— m — 1). Moreover, fc can change by no more than 

• Case I: the vertices of sampled edge are in the same n(Ci + Cj)-k and f a (t,u, fc), so we have that in one step: 
component; 

• Case II: they are in two components C and C, and 

AA; = min(n(Cj + Cj)-fc,/ a ), (2) 

Cj + Cj < fc/n; 

A* = -Am = <5(n(Cj + Cj)-fc,/ a ). (3) 

• Case III: they are in two components Cj and Cj and 

C + Cj > k/n. Here <5(x, y) = 1 if a; < y and 5 (or, y) = otherwise. 

According to the BFW(a) algorithm, we sample a ran- As u increases, the evolution equations of fc, t and m 
dom edge at step u: in Case I, the edge is also accepted; are established as follows: 
in Case II, the edge is also accepted and two components 

Ci and Cj merge together; in Case III, we should consider dfc v-^ , , n n , 

— = 2 2^ min \ n \ C i + C j)- k > fa) CiCj, (4) 
the constraint condition t/u < a + l/\/2fc (on the 14th Ci+Cj>k/n 

dm 

line of the BFW(a) algorithm). This constraint condition = ~ F2 \ t > M ' k > 

is the kernel hard core of BFW(a) model, which ensures 

-2 ]T SMd + CA-kJJCiCj, (5) 

Ci + Cj>k/n 

either components are increasing evenly or dramatically. d£ 

— =P 1 (t,u,k) + P 2 (t,u,k) 
au 

Let's first introduce a function f a (t,u,k), which dc- 

+2 ]T 5{n{C i +Cj)-k,f a )C i Cj. (6) 
notes the maximum acceptable value of Afc at one step. d+Cj>k/n 

Due to the BFW(a) algorithm, if the rate of accepted Here the function Pi(t,u,k) is defined as the probability 
edges t/u is smaller than a, any sampled edge should that the vertices of a randomly sampled edge at step u 
be accepted; else, k can only increase until the condition are in the same component (Case I); similarly, the function 
t/u < a + \/\p2k is invalid. Thus, the function f a (t, u, fc) P2{t, u, k) is defined as the probability that they are in two 
is shown as follows: components with sum smaller than fc (Case II). Therefore, 

we can simply obtain: 



fa(t,u,k) 



min{* eN +|A> a+ - 7 =i= } ,ifA> a 

m 

if£<a. P 1 (t,u,k) = J2cf, (7) 



oo, 

(1) 



i=l 



P 2 {t,u,k) = 2 C i C r ( 8 ) 

According to the definition of f a , when a randomly Ci+Cj<k/n 

edge is sampled between two components Ci and Cj, if For Eq.(4), k is the upper bound of size of the largest 

and only if f a > n(d + Cj) — fc, we can accept the edge component and never changes in Case I and II; only in 
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Case III, k can change by no more than n{Ci +Cj) — k and For Pi (t, u',k) > a with some u', we can prove that 

f a (t, u, k), i.e., Ak = min (n((7j + Cj) — k, f a (t, u, k)). For t/u > a + l/y/2k for any u > u' . Doing calculations on 

Eq.(5), Am = in Case I and —1 in Case II respec- both sides of Eq. (6) from an initial state (to,«o) to a cur- 

tively; in Case III, the number of components will de- rent state (t\,ui), we have: 
crease by 1 if and only if f a > n(d + Cj) — k, i.e., 



-to = [ 

J u 



Am = —5 (n(Ci + Cj) - k,f a ). For Eq.(6), the sampled 
edge must be accepted in Case I and II, so At = 1; similar 

+ / 2 6 ( n ( C i + C o)~ fc > fa)CiCj du 



til 

Pidu 

Uo 

111 



to the m of Eq.(5) in Case III, edge can be accepted when Ju c +c >k/n 

f a >n(C i + C j )-ka,ndAt = 5(n(C i + C j )-k,f a ). > / p 1 du. 

J Uo 

No matter how large u is, we can always find some 
2.2 Steady-state conditions of evolution of the u > u with t a /u > a because BFW(a) model ensures 

BFW(a) model the fraction of accepted edges never smaller than a. Sup- 

posed t/u is always smaller than a, by the rule of BFW(a) 



Although the Eq.(4)-(6) are unsolvable, some interest prop 
erties and results can still be deduced from these equa- 
tions, especially the steady-state conditions. 



model, the sampled edge must be accepted and t, u in- 
crease accordingly in each step, which will lead to the 

increase of the value of t/u and finally make t/u > a. 
Taking the right side of Eq.(4) and (5) as 0, we can For ^ Qmy when ft sampled edge Hnking c , and c . 

obtain P 2 (t, u, k) = 0, min (n(C\ + Cj) - k, f a ) = and ig acceptedj the part of C 2 + q changes to {Q + Cj) 2^ 

S (n(Ci +Cj)-kJ a ) = 0. Furthermore, for two compo- which wiU make the yalue of p ^ ^ ^ fc) increase; otherwisej 

nents with sum smaller than k, they must merge together, p ^ ^ ^ fc) ^ neyer change _ go once Pi ^ u ^ k) > a for 

but merging operation is forbidden after system stabilizes, gome % it wiU bg kgpt for eyer 

so any two giant components stay with n^ + C,)- k>0 In summarVj once ^ k) > a for some „/ } choosing 

after stabilizing, then the steady-state conditions can be Uq > u > ^ ^ > ^ wQ obtain; 
simplified to be: 

t\ > t + a / du > au\. 

f a (t, u, k) = ■'"o 

(9) 

p 2 (t u k) = Notice that the formula above is correct for any u\ > u$. 

After the giant components come up, we have: 

According to the definition, f a (t,u, k) = if and only 



if t/u > a + 1/V2fc, so we just need to prove t/u > a + 

\/\p2k and Pi(t,u,k) > a are equivalent. k <~ nC max ~ 0(n). 
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So l/\2k ~ o(l) and it can be ignored when n is large 
enough. Therefore, when P\(t, u\,k) > a for any u% > u', 
we obtain ti/ui > a + 1/ V2k. 

In the other side, if f a (t, u, k) = keeps for any u > u' , 
Eq.(6) turns to be: 



dt 



— = P 1 (t,u,k). 
du 



In order to ensure t/u > a + o(1)(Vm > u'), we need the 
slope Pi(t, u, k) > a.D 

Finally, we can obtain the steady-state conditions for 
any u > u': 



(10) 



Pi(t,u, k) > a 
P 2 (t,u,k) = 0. 
In fig.l, simulations have verified this conclusion. For dif- 
ferent values of a, the values oit/u and k/n change until 
Pi{t,u,k) = £™°<7? > a > which means the BFW(a) 
system evolves until Eq.(10) taking effect. 

2.3 Merging mechanism on the multiple giant 
components of BFW(a) model 

Since in Case I and II, sampled edges are always accepted; 
we just need to explicitly consider the Case III. Consider- 
ing t/u ^ a + 1/V2k at step u; according to the BFW(q) 
algorithm, when the edge e n +i is sampled at step u + 1 
and t/(u + 1) ^ a + l/v / 2fc, the edge e u +i is simply re- 
jected. But if t/(u+ 1) < a + 1/V2k, k needs to increase. 
To obtain the maximum accepted change of k, i.e. f a , we 
set t/u to be the smallest value: 





a=1/3 


t/u 






g(k) 






P,(t,u,k) 






k/n 








Fig. 1. Simulation of BFW(a) model by 100 random in- 
stances with N — 10 6 nodes. The solid lines denote the frac- 
tion of accepted edges for different values of a, which can 
be calculated by t/u in algorithm; the dashed lines denote 
Pi(t, u, k) = X/i=i Cf- The above subgraph shows the dynami- 
cal evolution of variables k/n, t/u and Pi(t, u, k) with a = 1/3. 
By the results, when Pi(t, u, k) > a, t/u keeps greater than a 
and k keeps still. 



2k 



differentiating u on both sides by k we find that: 



du _ 1 t 

dk ~ ~ k( \ /r^rr\ 2 k^ 



2y/2 



(11) 



Before the giant component appears, we consider that 
t ~ 0(n), k ~ nd ~ o(n). With du = 1 and Eq.(ll), we 
obtain f a ^dk^ 0(ki /n). Let S denote the component 
set {Ci I k/2 < Cin < k}. If the sampled edge links com- 
ponents Ci and Cj, where Ci, Cj € S, then n(Ci + Cj) — 
k - 0(k) > 0(0 /n). That leads to 8 {n(d + C 3 ) - k, f a ) 
0, which means this edge is rejected. So the edge linking 
two components of S must be rejected and only the edge 
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Fig. 2. An example for the evolution of BFW(a) model. The 
blue (yellow) components are (not) in the set S. The green 
cycles denote the isolated nodes. Subgraph (a) provides the 
initial graph, in which two random edges (ei,e2) are picked in 
each step yet only one is added to the graph based on selec- 
tion rule of BFW(q) model, whereas the other is discarded. 
Subgraph (b) is the situation of adding edge ei, two tiny com- 
ponents merge into a giant one. Subgraph (c) is the situation 
of adding edge e2, a tiny component merges into a giant one, 
forming a new member of set S. In subgraph (d), the sampled 
edge ez linking two giant components is rejected. 

linking to at least one component of CiS can be received. 
For the process is performed successively, either a new 
member in S comes up or the scale of an original one in 
S becomes more close to k. That is the key for coexisting 
multiple giant components, and they are expected to grow 
simultaneously before a critical point(fig.2). 

Similarly, after the giant components appear, the k ~ 
nCi ~ 0(n),VCi S S, so we have: 



f a (t,u, k) ~ 0(n 1/2 ) < n(d +Cj)- k. 



(12) 



Then we still have 5 (n(C + Cj) — k, f a ) = and the 
members of C<S keep merging into S until P2 = 0. Ac- 
cording to the steady-state conditions Eq.(10), if C sat- 
isfies Pi > a, the system will stabilize; else, Pi(t,u, k), 
which is the probability that the vertices of sampled edge 
are in the same component, is smaller than a. As it is 
proved above, if Pi < a, we can always have some u with 
t/u < a + 1/ v2fc, which makes k keep increasing (the 
14th and 15th lines of the algorithm) until two compo- 
nents merge together. 

Furthermore, only two minimum components (marked 
as C{ nm ,C™ m ) can merge together. We define P as the 
probability of any other two components (marked as Cj , Cj ) 
merging together before the two minimum ones, then 



P ^ (l - 2C7 m C^ m ) 



yn(C i +C i -C? in -C? in )/ f a (t,u,k) 



(.13) 



For one step, Afc = min (n(Ci + Cj) — k,f a ) = f a . Let's 



take f a (t, u, k) as average increase of k for one step. Based 



on Eq.(12), f a (t,u, k) ~ 0(n 1 ^ 2 ). When k increases to be 
larger than n(C^ m + C^ m ) but smaller than n{C l + Cj), 
that only one edge linking C™ m and C™" is sampled can 
make them merge together. So in order to ensure that 
Qmm an( ^ (jmm can '^ me rge together before Ci and Cj, 
we need k to increase by n(Ci+Cj — C™ m — C™ m ) without 
any components merging, which means n(Ci+Cj — CY lm — 



CT n ) I ' f a {t, u, k) edges not linking Cj™" and C^ in should 
be added. So the probability P satisfies Eq.(13). 
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Fig. 3. Up: The maximum component's size in critical interval 
of the BFW(a) process; the horizontal axis t reflects a "jump 
point" at which the maximum component increases dramati- 
cally. Down: the yellow and red rectangles represent the first 
and second minimum components before "jump point", and 
the blue ones represent the maximum component after "jump 
point". At the "jump point", the two minimum components 
emerge into the maximum one. All the other components in S 
keep unchanged at "jump point". 

As the system size n — > oo, the number of needed edges 



n(Ci + Cj - Cr™ - C% lin )/f a (t,u,k) -> 0(n x / 2 ) -> oo, 
causes P — > 0, which means two minimum components 
can merge before any other two components. This phe- 
nomenon can also be verified by simulation of BFW(cv) 
model (fig. 3). 

3 Quantitative properties on the giant 
components of BFW(a) model 

Since we have analyzed the BFW(a) model with fixed a 
in detail, in this section, we are going to calculate the 
number and size of giant components with arbitrary a in 
theory. Taking C™(cv) as the ith largest giant component 
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Fig. 4. Size of giant components in set S with different a. 
cycles, triangles, squares, diamonds, stars and crosses denote 
the fraction of giant components in steady state of m = 2, 3..., 7 
respectively. The results are obtained by 100 random instances 
with fO 6 nodes. The above subgraph denotes the scaling of the 
error of all | C 1 m+1 - CYU \ for every m = 2, 3..., 7 and i G 
[1, m — 1]. Here the C™ is average size of ith largest component 
in m-steady-state by simulations. 

of steady-state with m giant components (m-steady-state) 
to replace Ci above, we will find some common properties 
of C^ipc) when a belongs to some intervals. 

3.1 Steady states with different evolution parameter a 

fn BFW(a) model, parameter a plays a key role on the 
problems when the system can stabilize and which state 
the system can stabilize in. Defining (a m +i,a m ] as the 
m-steady-state interval in which there exists m giant com- 
ponents, 77i == 1, 2, 3 • ••, for any a with cv m +i > a ^ ev m , 
the value of Cj m (a) are all the same when system stabi- 
lizes (fig. 4), so we take C™ instead of C™(a) briefly for 
all a S (a m +i,a m ]. Thus, the members in set S evolve 
similarly with different phases in different intervals of a. 
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Moreover, if system BFW(a) with a £ (a m+ i, a m ] has 
m+1 components CJ™ +1 (a) in S, there must be Pi (t, u, k) < 
a, which leads to its collapse and a steady phase of m- 
steady- state. As mentioned above, Eq.(13) ensures that 
only the two minimum components can merge before sys- 
tem stabilize and they merge to the largest one in m- 
steady-state (fig. 3). So the C™ and C™ +1 must satisfy: 




Fig. 5. Critical points of phase transition of stable giant com- 



C£ 1 =Cr +A ,Vi = l,2,...,m-l. 



(14) 



ponents' number for a £ [0.05, 1]. The red solid cycles and 
In addition, when two components merge together, black lines represent numerical results of 100 random instances 



P 1 (t,u, k) will "jump" by (C™+1 + C™+ 1 ) 2 - (C™^ 1 ) 2 



with 10 6 nodes. The blue lines denote theoretical results. The 
solid triangles in below subgraph denote the error between nu- 



(C™ +1 ) 2 = 2C™ +1 C*™j 1 1 . Notice that in m- steady- state, 

. . ii-i / uicrical and theoretical results. 

Pi(t,u,k) keeps unchanged and is larger than a £ {a m +i , &m\ > 

so Pi(t, u, k) must be the upper bound of a in m-steady- i s expected to be the sampled edge number u. As rejected 



state: 



(15) 



Suppose all components' size in S are very close and 
2<=i Cj — 1) we have theoretical expression of the a m : 



= V(cd 2 = = -• 



1=1 i=l 

As the value of a goes smaller, the assumption is more 
close to the truth by numerical results (fig. 5). 



edges are almost between two components of S, the size 
of whole set S is almost unchanged if we take these re- 
jected edges on. According to the method of generating 
function |33) . the fraction x of giant component satisfies 
the equation 1 — x = e - 2xu / n ; here u represents the num- 
ber of added edges. At the critical point, the threshold 
t c ~ 1 and t c /u ~ a when u is large enough (fig.l). With 
(16) Eq.(15), we obtain the general equations of the fraction of 
giant components (7™ for any integer m: 



E™i(CT) 2 = « m 

YZiCT = xm ( 17 ) 

1 _ T — ~2x m /a 

Firstly, we can take the whole set S as a component and Here x m denotes the fraction of whole set S which has m 



3.2 Number and sizes of multiple giant components 



assume it grows similarly as the giant component on Erdos- giant components. Notice for case m — 2, Eq.(17) can be 
Renyi random graph, in which the number of added edges solved uniquely (due to error of a m in fig. 5, the accurate 
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result a m = 0.52 is adopted), then with the results and 4 Conclusion and Discussion 



Eq.(14), we can obtain all the sizes of multiple giant com- 
ponents. Contrast between theoretical and simulation is 
showed as follows: 

Table 1. The fraction of components, Pi and fraction of the 
whole set S with a = 1/2,1/3,1/4,1/5. The results are ob- 
tained by 100 random instances with 10 6 nodes. The num- 
bers in brackets are theoretical results, which are obtained by 
Eq.(14)-(17). 



We detect the steady-state and evolution process of BFW(a) 
model with both numerical and theoretical methods. Ac- 
cording to the rule of BFW(a) model, function f a is de- 
fined to calculate the change of the stage k and num- 
ber of giant components m. Furthermore, by establish- 
ing the mathematical expression of evolution equations on 
this model, an equal relationship between the parameter 
a and steady-state condition is proved. Meanwhile, with 
some hypothesis, the correspondence between parameter 
a and the number of giant components in steady-state 
is obtained, that is when a € (^qrf, BFW(a) model 
must stabilize with m giant components. Through the fur- 
ther analysis of the evolution process and the numerical 
results, set S is defined to find the rule of two components 
merging before and after the threshold. Moreover, sizes 
of giant components for different evolution strategies also 
has a close connection with each others and satisfy some 
constraint equations, which is derived from the evolution 
equations. 

So far, we can calculate the number and sizes of giant 
In summary, parameter a determines when the system components for different evolution strategies with theoret- 
can stabilize and which state the system can stabilize in. In ical methods, which can correspond with simulations very 
the evolution process of BFW(a) model, a can only take closely, especially when the value of a is smaller than 0.25. 
effect on when to increase k. As to how much k increases, a Additionally, the analysis of the steady-state and evolu- 
doesn't work. With this special evolution rule of BFW(a) tion process is of great help to explain why the percola- 
model, the connection between two adjacent steady-states tion of BFW(a) model is explosive and how explosive it is, 
is found and sizes of giant components are obtained. As which are almost supported by simulations before. For ex- 
the value of a decreasing, theoretical results can be much ample, before the the giant component appears, we can ob- 





a = 1/2 


a = 1/3 


a = 1/4 


a = 1/5 


Ci 


0.5736 


0.4142 


0.3220 


0.2631 


c 2 


0.4144 


0.3217 


0.2629 


0.2223 


c 3 




0.2633 


0.2223 


0.1926 


c 4 






0.1928 


0.1699 










0.1519 


EC 2 


0.5007 


0.3444 


0.2594 


0.2077 




(0.5000) 


(0.3333) 


(0.2500) 


(0.2000) 




0.9880 


0.9992 


0.9999 


0.9998 




(0.9802) 


(0.9975) 


(0.9997) 


(0.9999) 



better verified by simulations. 



taindfc/di ~ 0(k^ /n), although the k ~ 0(nz +s ) <~ o(n), 
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for < S < 1/3, we still have dfc/di ~ oo; here k ~ C max , 8. R. M. Ziff, Phys. Rev. Lett. 103, (2009) 045701. 

so the result means percolation of BFW(a) model is ex- 9. R. M. Ziff, Phys. Rev. E. 82, (2010) 051105. 

plosive. 10. O. Riordan and L. Warnke, Science 333, (2011) 322. 

Besides, we just analyze this model on Erdds-Renyi 11. R. A. da Costa, S. N. Dorogovtsev, A. V. Goltsev, and J. 
random graph; as to other random graph with any degree F. F. Mendes, Phys. Rev. Lett. 105, (2010) 255701. 

distribution, the theoretical methods in this paper also 12. P. Grassberger, C. Christensen, G. Bizhani, S.-W. Son, and 
work, which only need to modify the related probabilities M. Paczuski, Phys. Rev. Lett. 106, (2011) 225701. 

of the vertices of randomly sampled edge are (not) in the 13- H. K. Lee, B. J. Kim, and H. Park, Phys. Rev. E. 84, 
same component. (2011) 020101(R). 

14. L. Tian and D.-N. Shi, Phys. Lett. A 376, (2012) 286. 

mi ii . 15. M. x. Liu, J. f. Fan, L. s. Li and X. s. Chen, Eur. Phys. J. 
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